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(^ [ Abstract 

The method of constructing the tomographic probability distributions 
Mh, describing quantum states in parallel with density operators is presented. 

Known examples of Husimi-Kano quasi-distribution and photon num- 
C I ber tomography are reconsidered in the new setting. New tomographic 

schemes based on coherent states and nonlinear coherent states of de- 
formed oscillators, including 5— oscillators, are suggested. The associated 
identity decompositions providing Gram- Schmidt operators are explicitly 
given. 

Key words Quantum tomograms. Coherent state tomograms, Photon num- 
ber tomograms, Husimi-Kano quasi-distribution, Sudarshan's diagonal co- 
j^ ' herent state representation, f— and g— oscillators. 

PACS: 03.65-w, 03.65. Wj 

1 Introduction 

There are several representations of quantum states providing the possibility to 
present equivalent, but different in their form, formulations of quantum mechan- 
ics pp. The representations of quantum states are based on different integral 
transforms of the density operator [51 [3] taken in the position representation. 
The density operator in the position representation is mapped by means of the 
integral transforms either to Wigner quasi-distribution function [4j , or Husimi- 
Kano iiT— function [3 [6]. In this paper we have decided to keep up with the 
original notations of the pioneer papers on the subject. Quasi-distributions are 



usually referred to as phase space representations of quantum states. Another 
important phase space representation is related to the Sudarshan's diagonal co- 
herent state representation [71 [5] . Recently the tomographic representations of 
quantum states was suggested UHl E] using the Radon integral transform of 
their Wigner functions. The tomographic representation exhibits some specific 
property in comparison with the other phase space representations. The tomo- 
graphic probability distributions (tomograms) associated with quantum states 
are standard positive probability distributions. The mathematical mechanism 
of constructing the tomographic probabilities in abstract Hilbert spaces was elu- 
cidated in [l^ for the finite dimensional case. Also infinite dimensional Hilbert 
spaces were considered in the context of mapping vectors and density operators 
onto tomographic probabilities in [13] . 

The details of those tomographic approach constructions inspire to study 
other possibilities to know better and create new schemes of description of 
quantum states by positive functions. First of all, the Husimi-Kano ii'— function 
interpreted as quasi-distribution in the phase space representation is a real non- 
negative bounded function of two variables, q and p. Strictly speaking, these 
variables cannot be interpreted as the actual position and momentum since the 
uncertainty relations prohibit the existence of a distribution in the phase space 
depending on variables which cannot be measured simultaneously. Nevertheless, 
being a real non-negative and normalized function on "phase space" , it has a 
chance to be interpreted from the point of view of the tomographic approach. 
One of the aims of the present work is to extend our analysis of the mathematical 
mechanism of constructing tomographic sets to extra situations, by including 
into the construction properties of the coherent states [14] and properties of 
i^— function and Sudarshan's diagonal coherent state representation [15] in the 
sense of their relation to distributions (i.e., generalized functions) . We call this 
extension coherent state tomography. From this point of view the Husimi-Kano 
iiT— function will be interpreted as a tomogram of a quantum state. 

The coherent states are closely connected with linear vibrations (linear har- 
monic oscillator). The deformed oscillators, e.g. q— oscillators [TBlfTT] and /— 
oscillators [TS] dH] are related at a classical level with specific non-linear vibra- 
tions, so that non-linear coherent states were introduced [501 HI] to describe 
the corresponding states of a non-linear quantum oscillator, which yields in the 
linearity limit the standard coherent state counterpart. Another goal of this 
work is to develop the /—deformed version of the tomographic approach and 
in particular to present, for a specific choice of the non-linearity coded by the 
function f = fq, the g— deformed coherent state tomography. 

Also, by using /—deformed coherent states, we develop the /—deformed 
version of the photon number tomography [551 HSl HI] as well as we obtain 
the identity decomposition for the deformed tomographies. In fact, we are 
addressing the general problem of constructing a map from a Hilbert space (or 
projective Hilbert space) onto a family of probability distributions (tomograms). 

As in the usual formulation of quantum mechanics there are several schemes 
like Schrodinger picture, Heisenberg picture, Dirac picture, even in the tomo- 
graphic approach exist several different schemes, like symplectic tomography, 



photon number tomography, optical tomography and so on. Our last aim is to 
point out the common general mechanism of constructing all these tomographic 
schemes and in this way to extend the list of the tomographies by including the 
deformed ones. 

The paper is organized as follows. In section 2 a short account of our con- 
struction of a tomographic setting in abstract Hilbert spaces is given. In sec- 
tion 3 coherent state tomography is developed. In section 4, a resolution of 
the unity is obtained which provides a relation between Sudarshan's diagonal 
coherent state representation and Husimi-Kano iiT— function in the suggested 
tomographic setting. That identity resolution is then recovered in the context 
of the generalized phase space distributions associated with the Agarwal-Wolf 
ri— operator ordering [25l [26]. In section 5, generic deformations of coherent 
states connected with /—oscillators, g— oscillators, and s— deformations asso- 
ciated to operator ordering are developed. In section 6, photon number to- 
mography is reviewed and the new scheme of deformed tomography with both 
nonlinearity and s— ordering is extended to photon number states. Some con- 
clusions and perspectives are drawn in section 7. An Appendix presents a check 
of both the resolutions of the unity associated with the coherent state and the 
photon number tomographies and a derivation of the Mehler's formula in the 
holomorphic representation. 

2 Preliminaries 

In this section we give a short account of our previous works [121 113| , in which 
we have given an interpretation of quantum tomography in an abstract Hilbert 
space Ti by means of complete sets of rank-one projectors {P^i} ^^^ , where M 
is a set of (multi-) parameters, discrete or continuous, collectively denoted by /j,. 
In general, a tomogram of a quantum state {ip) is a positive real number 7^(/Lt), 
depending on the parameter /i which labels a set of states |/i) e H, defined as 

T4^,)■.^M4>)\\ (1) 

Our main idea was to regard the tomogram %p{iJ-) as a scalar product on the 
(Hilbert) space H of the rank-one projectors |^) (/i| = -Pp ^ |Pp) G H : 

T^{^l)^^r{P^p^)^■.{P^\p^). (2) 

Equation ^ may readily be used to define the tomogram of any density operator 
p or any other (bounded) operator A 

TA{fi):=Tr(p,A)={P,\A). (3) 



Equation ^ shows in general that to any operator A a function (/i|A|/i\ of 

the variables /i corresponds in a given functional space. So, a tomograph may 
be thought of as a de-quantization, and in fact is an useful tool to study the 



quantum-classical transition by comparing classical limits of quantum tomo- 
grams with the corresponding classical tomograms f33^. In the same sense, the 
inverse correspondence T^Cm) —* A may be considered to give a quantization. 
The reconstruction of the operator A from its tomogram Ta (/i) may be written 
as 

A=Y. ^mTi [p,a) ^\A)^Y. IG'm) {PM) • (4) 

In other words, the reconstruction of any operator is possible because the to- 
mographic set {P^} p ,, provides a resolution of the identity (super-) operator 

onH: 

We may then view jG^) and {P^\ as dual supervectors. For instance for the 
spin tomography, in the maximal qu-bit case M is the Bloch sphere S*^ of all 
rank-one projectors and for any operator A holds [12]: 

A = G{e,<l))TT{P{9,(l))A)sin9d9d(l> (6) 

Jo Jo 

27r pTT 

/ P{9,(f>)Tr{G{e,<j))A)sineded(f>, (7) 

Jo Jo 

where, in matrix form, 

, , _ 1 r l-hcose* e-*"*sin6l 1 /V/^ ,n _ J_ T 1 + ScosS SC't'sme 
^(fc'.W - 2 [ e*"*sin6l 1 - cosg J ' ^ ' "^^ ~ 47r [ 3e*"*sin6i l-3cos6' 

By definition, the set {Pfj,} ^j^^ is tom,ographic if it is complete in H. A 
tomographic set determines a tom,ograph which is a functional, linear in the 
second argument 

(P^,i)^Xi(M) = Tr(p,,i). (8) 



This definition is appropriate in the finite n— dimensional case, where 

l/i) enn^P^.e H„2 = BiHn) = Hn ® H„, 

but in the infinite dimensional case the relation EI — B{Ti.) is no more valid and 
there are several relevant spaces, as the space of bounded operators B{T-C) and 
that of compact operators C(7i), the space of Hilbert-Schmidt operators ^2 and 
that of trace-class operators 3i. Their mutual relations are: 

3i C J2 C C{n) C B{n). 

B{Ti) (and G(7i)) are Banach spaces, with the norm ||A|| = sup(||,|/,||^]^-) ||^'0I1 i 
while ^2 is a Hilbert space with scalar product {A\B) = Tr (^A^ B). Finally, 3i 
is a Banach space with the norm ||A||-^ = Tr(|A|) . The following inequalities 
hold true: 

ll^ll< 11^112 < Pill- 



So ^2, the only Hilbert space at our disposal to implement our definition of 
tomographic set, is endowed with a topology which, when restricted to the 
trace-class operators, is not equivalent to the topology of 3i. This may have 
serious consequences. In fact, in the finite dimensional case, the set {P^} ^^^ 
is complete iff 

Tr {Pf,A) = y^le M =^ A^O. (9) 

Such a condition guarantees the full reconstruction of any observable from its 
tomograms. Now, in 32, Eq- ® reads: 

{P^\A)^0 yfie M =^ A^O & AeJ2- (10) 

Then, as 3^2 is a *— ideal in B(Ti.), there may exists a non-zero operator B, which 
is bounded but not Hilbcrt-Schmidt, such that 

Tr {Pf,B) = V^i e M 

In other words, different observables may be tomographically separated only 
when their difference is Hilbert-Schmidt. Nevertheless there is a second case, 
when the set {Pfj,} ^^ of trace-class operators is complete even in 3i. Then, 
recalling that 3i is a *— ideal in its dual space B{T-i): 

31 = Bin), 

the expression Tr {Pf_iA) is nothing but the value of the linear functional Tr {-A) 
in P^. Hence, Eq.® holds unconditionally 

Tr(P^yl)=0 yfieM=^0 = \\Tr{-A)\\=\\A\\=^A = 0. (11) 

Thus, the finest tomographies are those based on sets of rank-one projectors 
which are complete both in ^2 and in 3i. As a matter of fact, this is the case 
for the main tomographic sets, like the symplectic, the photon number and the 
coherent state tomographic sets. 

From a geometrical point of view, tomographic sets are "skew " sets of 
projectors. In other words, we define a set of projectors to be skew when it 
spans the whole Hilbert space. Thus, any tomographic set is skew as it is 
complete. Besides, we define a set 7 of projectors to be locally skew in Pq , if it 
contains Pq and any neighborhood of Pq contains a skew subset of 7. Perhaps 
the simplest case of a tomographic set which is skew and locally skew is provided 
by the coherent state tomographic set discussed in the next sections. 



3 The coherent state tomography 

This tomographic set is generated by the displacement operators {V (z)} de- 
pending on a complex parameter z 

X'(z) = exp(za^-z*a) , z e C, (12) 



which acting on the projector |0) (0| of the vacuum Fock state, a|0) — 0, yield 
the projectors 

|z)(z|=I?(z)|0)(0|I?(^)^ ^eC, (13) 



associated to the usual coherent states 



\z) 



I |2 I |2 oo j 

exp(-^) exp (^at) exp (-z*a) |0) = exp(-^) ^ -^a^^ |0) 



We recall that the coherent states are a (over-) complete set in the Hilbert 
space Ji. Any bounded set containing a limit point zq in the complex z— plane 
defines a complete set of coherent states containing a limit point, the coherent 
state l^o) ) in the Hilbert space H.. In particular, any Cauchy sequence {zfe} of 
complex numbers defines a Cauchy sequence of coherent states {|zfe)} , which is a 
complete set. The same holds for any extracted subsequence. This completeness 
property holds as exp ( \z\ /2 j (z|V') is an entire analytic function of the complex 
variable z*, for any ji/;) G 7i. Then 

(zfe|V;)-0 Vfc^|^)-0, (14) 

because Zq is a non-isolated zero of an analytic function. 

Besides, any bounded operator A may be completely reconstructed from its 
diagonal matrix elements (z^ \A\ Zk) ■ In fact, exp ( \z\ /2 + |z'| /2) (z \A\ z') is 



an analytical function of the complex variables z* , z' , so it is uniquely determined 
by its value exp f \z\ j (z \A\ z) on the diagonal z' = z. This is an entire function 

of the real variables SRz, 3z, which is in turn uniquely determined by its values 
on any set with an accumulation point. 

The rank-one projectors associated to a complete set of coherent states are 
complete in the Hilbert space 32- In particular, any Cauchy sequence {|zfc)} 
generates a tomographic set {\zk) {zk\}- In fact, bearing in mind the previous 
remark on the reconstruction of a bounded operator, it results 

Tr(yl|zfe)(zfe|) = (zfe|A|zfe)=0 Vfc ^ A - & AeB{n). (15) 

This shows that a tomographic set of coherent state projectors is complete even 
in 3i. So it is skew. Moreover, any extracted subsequence (Iz^;^) (z^^l} is again 
complete, so {\zk) {zk\} is locally skew in its limit point. In the next section we 
address the case when z varies in the whole complex plane. 

4 The resolution of unity generated by the co- 
herent state tomographic set 

It is possible to interpret the well known Husimi-Kano iiT-synibol of a (bounded) 
operator A as the coherent state (CS) tomogram of A: 

Ka{z) :- (z |i| z) =: Tr(|z> (z| A). (16) 



In particular, when A is chosen as a density operator p, the identity holds 

-(z|p|z) = Tr(p) = l, (17) 



which allows for the probabilistic interpretation of the CS tomography. As 
a matter of fact [n| the /("-symbol exists also for a number of non-bounded 
operators. The CS tomographic set is complete both in ^2, the space of Hilbert- 
Schmidt operators, and in 3i, the space of trace class operators acting on the 
space of states. In fact, the formulae 



and 



A 



A 



,\.T 



(Pz (fz' 



A 



l^> {A 



(18) 



= e 
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E 

n.TTi— 



(z*Y{z'y 



d 



\n-\'m 



az*"(9z™ 



(el^l^(z|i|z)) 



show that if the tomograms 



A 



(19) 



z ) of a bounded operator A vanish for any 

z e C , then A is the zero operator. Equation p9p is implicit in Eq. (6) of 
Sudarshan Ref. [7]. 

So, a resolution of the unity exists, which allows for the full reconstruction 
of any (bounded) operator from its CS tomograms. We are interested in the ex- 
plicit determination of such a formula. Now, the Sudarshan's diagonal coherent 
state representation 4'j^{z) of an operator A is defined through the equation 



A= —cj,^{z)\z){z\. 



(20) 



Bearing in mind the previous qu-bit reconstruction formula, Eq.([7]), we may 
note the analogy with Eq. pp)) defining 0^(z). Then, if our guess is right, we 
need to get explicitly the first form of the reconstruction formula, Eq.®. In 
other words, we have to invert the well-known relation: 



Ka{z') 



A 



—<I^A{z)\{z\z')f 

TT 



(fz 



(t>Aiz)e 



(21) 



which follows at once from Eq. ([20|l defining 4>a{z)- This relation shows that 
Ka{z') is given by the convolution product of 0^ times a gaussian function. 
Then, denoting with Ka{z'j^, Zj) and 0^(zfl,z/) the K and symbols, with 
z' = z'j^ + izj : z = Z]^ + izj, the Fourier transform [27] of Eq. pT|) reads: 



(J/Z ry i-iiZ J 



2-K 



X^(z^,z;)e-^(«^«+''-i) 



(22) 



dz'fidzj f dzudzj 

27r 



-0^(z;?,z/)e-(^«-^«)'-(^^-^^)'e-^(«^«+''^^) 



and we readily obtain 



from which 






that formally yields 



dS^drj 



iew)/^ 



KA{tv)e 



i(izn+vzi} 



(23) 
(24) 

(25) 



The presence of the anti-gaussian factor shows that the inverse Fourier transform 
of 0^(C, Tj) exists only when the asymptotic decay of Ka{^, if) is faster than the 
growth of e'-^ +'' ^/^. However, the integral always exists as a distribution, as 
was proven in Ref. 15 . By virtue of this remark, we may go on and substitute 
the previous expression into Eq. pp)) getting 



A^ 



d^dr] 



27r 



■.^^'+^'y'KA{^,vy 



K^ZR + TlZi) 



(26) 



d^drj f dz'j^dzj 



27r 



27r 



ifA(4,4)e(«'+''')/V[«(^«-^«)+"(-^--0] 



Upon interchanging the order of integration, we may write the expected recon- 
struction formula as 

A^ j ^G{z')Ka{z'), (27) 

where the Gram-Schmidt operator G(z') reads; 

G(,') := / ^ / '^A^e^eW)/^M^.~^'n)M^'-^'.)] |.) (,| . (28) 

J 2tt J 2tt 

In other words, the resolution of the unity generated by the CS tomographic set 
is 



2^' 



d'z 



G{z')Tr{\z'){z'\-) 



(29) 



We observe that 

Kg{z-){z) = (z G{z') 



d^dr] 
( 

2tt 



,4-«^k-'^^;]e(«^+''^)/4i^|^^^^l(^,^) 



^^e'[-i^'.-vz',]l 



2tt 



^\z){z\{£„v) ^<l>\z){z\ (z') 



(30) 



This, substituted in the reconstruction formula Eq. ([77|) . gives also the repro- 
ducing kernel formula 



ifG(.")W- / -^{^ G{z') z) KGi.,,){z') 



-KG(z'){z)KGi..,){z') 



(31) 



or equivalently 



Since 



^kX-l (^") = / -;f'^l->(-l (^')0|.'>(.'l (^") • (32) 



1^) (^1 = / ^'^k)(.| (^') 1^') (^'1 ^ -^I.X.I (^') = '^'5(z - z') , (33) 
we remark that Eq. (|30|) amounts to the orthonormahty relations 

Xg(.')(^) = Tr(|z) (z| G(z')) = ^5(z - z') (34) 

between dual sets of supervectors. 

In view of Ea. (fT9|) . we can use the previous expression of Kq to check the 
reconstruction formula, Ea. (|27p . in matrix form as: 



(z \a\ z)^ j -^Kg(^,^{z)Ka{z') = Ka{z). (35) 

A check of the reconstruction formula in the position representation is given in 
the Appendix. 

In view of the generalization of our results in the next section, it is expedient 
to make preliminarily contact with the generalized phase space distributions 
associated with the Agarwal-Wolf fi— operator ordering [25, 26 . In fact, in the 
context of that operator ordering theory, the K and </> symbols, appearing in our 
direct and dual reconstruction formulae, are related to the Wick (i.e., normal) 
and anti-Wick ordering respectively. As it is well known, it is possible to express 
any operator A(a, a^) in terms of a function Fa{z, z*) of two complex variables 
u = z,v = z* , by means of an operator Cl which may be realized as integral 
operator 

A{a,a^)^ f—FA{z,z*)A^^\z-a,z*~a^) (36) 

where the Vl— ordered delta function 

A(")(z-a,z*-at)= / — fi(u;, u-*) exp(w;(z* - a^) - u;*(z - a)) (37) 

is obtained by the displacement operator and an entire analytic function n{w, w*) 
of two complex variables, which characterizes the particular rule of association 
between A and Fa- For instance, the Wick ordered delta function A^'^) is given 
by 0(w,w*) = exp(ww*/2), while the anti-Wick ordered A*^^-' is given by the 
reciprocal function exp(— wk;*/2). 

The function Fa is associated to the operator A via the anti-fi— ordered 
delta function: 

Fa{z, z*) ^ Ti-[A{a, a^)A^^\z -^ a, z* - a^)]. (38) 



It is apparent the similarity between Eq. ([55)1 and our previous reconstruc- 
tion formulae, but Agarwal-Wolf formulae, in general, do not contain rank-one 
projectors, which instead characterize our tomographic theory. However, it re- 
sults that the anti-Wick ordered delta function is a rank-one projector: 

A(^)(z-a,z*-a^)^|z)(z|, (39) 

so that the Wick ordered delta function A^^^ is nothing but our Gram-Schmidt 
operator G of Eq. ([28]) , and Eq. ((38l) yields the Husimi -Kano function Ka while 
Eq.(f36|) is just our direct reconstruction formula Eq. ([27|) . When exchanging 
the role of the two orderings one obtains respectively the Sudarshan's diagonal 
representation (j)^ and the dual form of our reconstruction formula, Eq. (|20p . 

5 The (/, s) —deformed CS tomographic set 

The usual coherent states, that we have recast in a tomographic setting in the 
previous sections, date back to Schrodinger. We have defined them by applying 
the unitary family of the displacement operators {I? (z)} on the vacuum state 
|0) , such that a |0) — 0. Recently, a number of so-called generalized coherent 
states (GCS) have been introduced, in connection with oscillator algebras dif- 
ferent from the standard one. In particular, an interesting family of such GCS, 
related with deformed oscillator algebras, is the set of nonlinear CS [IHIIIIIIIH] 
which arises from a deformation ruled by a class {/(«)} of functions of the 
number operator. This class is the set of real functions such that 

0<c<[/(n)]!<c'<oo, neN; [/ (n)]! := / (n) / (n - 1) .../ (1) . (40) 

Following '28], that we briefly recall mainly to establish notation, we assume 
that these functions satisfy some further condition: they are continuous and 
moreover 

/(n) = l, ,, = 0,-1,-2,... . (41) 

Now, introducing the selfadjoint /—deformation operator 

OQ 

Ef:^f{n)\{f{h-m)=:[f{h)]\ (42) 

TTl — 1 

we have 

£j^a£f = af (n) =: A; £j^a)£f ^ -j^a^ =: A\ . (43) 

Then, the /—deformed displacement operator is 

Vf (z) = £J^V (z) £f = exp fzi^ - z* A\ . (44) 

Note that A\ appears instead of A^ = f (n) a' in the deformed displacement 
operator. This A\ obeys the commutation rule ^,^41 = [a,a^] = 1- 



10 



Now, it is possible to compose the /—deformation with the s— deformation 
of operator ordering [29l [30], so defining the (/, s) —deformed displacement op- 
erators as: 



I5/.S (z) :— exp ( —szz* j Vf (z) — exp ( —szz* j exp f zAl — z*A] , (45) 



where z G C, — 1 < s < 1. 

The (/, s) —deformed displacement operators, acting on the vacuum, gener- 
ate the (/, s) —deformed coherent states: 

\z; /, s) := e(i+^)l^l'/27V,,^I?/ (z) |0) = e^'+'^^'^'/^N^jSj'V (z) £f |0) 

= e^i+^^l^l'/^TV.jf/ip (z) / (0) |0) = e(i+")l"l'/2iv,_^£:-i \z) (46) 

where the normalizing factor is 



N. 



^J 



E 



ki 



2fc 



ofc!([/(fc)]!r 



-1/2 



(47) 



Now, we are ready for the tomographic setting of such deformed states. We 
define the (/, s) —deformed ii'— function of an operator B as: 



Kf/iz):=(z-J,s 



B 



z;f,s 



(48) 



(z |e(i+^)l^l'/2iV,,/£7iBe(i+^)l^l'/2iV,,/£7i| z) =: 6^'+'^^^^" NljKsif) (z) 



where 



B{f):^Ej'BEj 



(49) 



So, the (/, s) —deformed iiT— function of an operator B is the ilT— function of the 
/—deformed operator B (/) times a known c— number. 

Thus, from the reconstruction formula of the non-deformed CS tomography: 



B{f)= —G{z)Kb,j){z), 



n 



(50) 



we get the reconstruction formula for the (/, s) —deformed CS tomography 
^—£fG{z)£jK^^f){z) = /^ [£,G(z)f,e-(i+^)l^l^7V-| 

^2 



B 



J^Gf-{z)K{^^{z). 



(51) 



Remembering the expression of G, Eq. ([28]) , we have explicitly the (/, s) —deformed 
Gram-Schmidt operator G^''^ as: 



Gf'^iz) = - 



-{l+s)\z\' 



cPz' 

~2^ 



a z 

"2^ 



d£,dr] 
( 

27r 



,ieW)/4^^[iiz'n-zn)+v{z',~z,)]^^ 1^,^ ^^,| ^^ 



dS^drj 



keW)I^^Wn-nYr,{-i--i)\£2 1^/. j^ ^^ l^^i. f^ ^1 g 
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We recall [55] that £f is selfadjoint and moreover SfSl = £"? is a bounded 
positive operator (with bounded inverse) which satisfies the intertwining condi- 
tion: 

A^£] = £]A] . (52) 

In other words, the operator Al is nothing but the adjoint of the /—deformed 
operator A with respect to a new scalar product in our carrier Hilbert space 7i: 

(^, V)/ := {v, £]i^) ^ {Aip, ^) ^ - (if, A^£]i>) = (^, i}v) ^ (53) 

So, a better notation would be A^f instead of A\ , however we shall use the 
latter one. 

The presence, in the expression of the deformed Gram-Schmidt operator, of 
the deformed projector 

£]\z'-J,s){z';f,s\£] (54) 

is reminiscent of the new deformed scalar product, and allows the representation: 



w;/,s 



Gf^^z)\w'-f,s) 



-2(1-1-8)1^1^ 



(55) 



d^dr] 
27r J 2-K 

When w' — w, Mve get 



(eW)/^^Wn-^n)+r,{^':-^.)] (^. J^ ,|^^ J^ ,) (^^ f^ ,|^'. f^ ,) 



K: 



\w;f,s){w;f,s\ 



(z') := (u;; /, s|z'; /, s)j (z'; /, s|u;; /, s) ^. = -K'i2|^./,,)(,^.^,,|£2 (^') 

(56) 



and 



K, 



f.s 



(w) 



^-2{l+s)\z\^ 



w;f,s 
d^dr/ 



Gf^^iz) 



■w,f,s 



2tt ^£j\wj,s){wj,s\ej y^'^) 



^ej\w;f,s){w;f,s\£j ^^^ ' 



(57) 



which may be compared with the analogous formula of the standard CS case, 
Eq.®. 

Let us consider now an important case of particular deformations, the q— 
deformations. They are defined by considering the so called q— oscillators (THl 
I17j . For q— oscillators, the analog of annihilation and creation operators a, a' 
are the operators dq and aj which obey the relation: 



dgdl - qd\dq = q"' 



h = d'd. 



(58) 



12 



where g is a c— number. When q -^ 1, the relation of Eg. ([55)) becomes the 
standard boson commutation relation. One can introduce q— coherent states 
studying the eigenvalue problem 

aq\z;q) = z\z;q) . (59) 

The properties of g— oscillators and g— coherent states were intensively discussed 
in the literature. As it was shown, e.g., in Ref. [21]. the q— deformations are 
a particular case of /—deformations, with a specific function fq coding the 
g— nonlinearity of the g— oscillator vibrations, which reads 



sinh An 



/9(^):-\/^^, g:-eV (60) 

One has: 

aq=afq{n). (61) 

The operator Uq of Eq. ([CT|) and its adjoint aj satisfy the relation of Ea. ([551) . As a 
result, all the previous formulae for the deformed coherent states are valid even 
for g— deformations. In particular, this holds for the reconstruction formula 
Eq. (|5ip and the associated resolution of the unity, in which we have only to 
replace f by fq. 

6 Photon number tomographies 

6.1 The photon number tomographic set 

The photon number tomography is generated by the irreducible family {a^a, {"D (z)}}, 
where the displacement operators {V {z)} act on the number operator a) a of 
the harmonic oscillator, yielding the family of selfadjoint operators T{z) : 

f{z)^V{z)a^aD{z)\ zeC (62) 

which has the spectrum of cl^ci, eigenvectors \nz) = 'D{z)\n) and associated 
projectors \nz) {nz\. 

We recall that the photon number projectors' set, containing the complete 
set of the coherent state projectors, is a tomographic set complete both in 0^2 
and 2^1. For the same reason, any Cauchy sequence {|nz/j) {nzk\} is locally skew 
in its limit point [13] . 

The whole set of photon number projectors generates the resolution of the 
unity: 

i = ^ / ^GA(n,z)Tr(|nz)(nz|.) . (63) 

n— 

The Gram-Schmidt operator Gx is given by 

A + l\" , . /A-P'*' 



G,(„,.) = _^ _ virM—^] i>'W. («) 
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Here A is a real parameter, — 1 < A < 1, which labels the family of equivalent 
kernels G\ (n, z) . 

The check of the matrix form of the resolution of the unity, Eq. (|S5)) , in the 
position representation is given in the Appendix. 

6.2 The (/, s) —deformed photon number tomographic set 

We can rephrase the section on the (/, s) —deformed coherent states, by letting 
the (/, s) —deformed displacement operators act on the eigenstates \n) of the 
harmonic oscillator, so generating the (/, s) —deformed photon number states: 

|nz;/,s) : ^ c^^+'^^'^'/^N,jVf (z) \n) = e^^+'^^'^'/^N.jSj^V (z) £f \n) (65) 
= e(^+^)\^\'/^N,j£j'V (z) [/ (n)]! \n) = e^'+'^^'^'/^N^j [f {n)]\£j' \nz) 

where the normalizing factor Nzj is again given by Ea. (|47p . 

Now, we define the (/, s) —deformed X— function of an operator B as: 



^/ '^^/ ' 



Kg (nz) :— (nz]f,s B nz;f,s) (66) 

= e(i+^)l^l'iV2^[/2(n)]!(nz 

=: e(i+«)l^l=iV|,,[/2(n)]!if^(,)(nz). 

where 

B{f):=£j'B£j\ (67) 

Again, the (/, s) —deformed isT- function of an operator B is the if— function 
of the /—deformed operator B (/) times a known c— number. Then, the recon- 
struction formula for any bounded operator B in the (/, s) —deformed photon 
number tomography reads: 

B = Y. -^£fGx{nz)£fKBU){nz) =: ^ / -^G{\nz)Kf/{nz). (68) 

n=0 '' n=0 •' 

with 

'^^''^"^^ " [P{n)]\N^ g/GA(n^)g/ 



e 



2 

2, J 

{i+s)\A' 4 /A + l\" ,„ , , /A-l\^/^ 



mi^pH^r^'^^^ 



[P{n)]\Nlfl-X^ \X-l, 
which is readily obtained by using 

£fV {z) a^aV^ (z) £f = £JVf (z) a)aDf (-z) = £}Vf (z) A\AV-^ (z) . (69) 
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7 Conclusions 

To conclude, we recapitulate the new results of this work. We have constructed 
for infinite dimensional Hilbert spaces the tomographic sets based on coherent 
states and Fock states, respectively, and proved decomposition of the identity 
operator acting on B{TC), the space of bounded operators on the underlying 
Hilbert space 7i, in terms of rank-one projectors of the tomographic set and 
associated Gram-Schmidt operators. 

The properties of linear and non-linear /—oscillators were used to develop 
new schemes of the coherent state tomography. The /—deformed tomographic 
projectors were shown to yield a resolution of the (super-) unity simply by 
means of a linear non-canonical operator which may be used to deform the 
scalar product definition of Ti. The particular case of the deformed oscillators, 
namely g— deformed oscillators, were studied as example of the coherent state 
tomography. 

We addressed the general problem of mapping an abstract Hilbert space on 
the set of fair probability distributions (tomograms) describing quantum states 
in the framework of probability picture of quantum mechanics. 

We hope to extend the obtained results to the case of multimode quantum 
systems and entangled states in future papers. 
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8 Appendix 

8.1 Checking the resolution of the unity for the coherent 
state tomography 

In the position representation the reconstruction formula Eq. ()27l) reads: 

(fz' 



A 



n 
dydy' 



G{z') x')Ka(z') 
£z' 



(70) 



G{z') x'){z'\v){v'\z') 



A 



y 



So, we check the resolution of the unity generated by the CS tomographic set 
by showing that 

<fz' 



G(z') x'){y'\z') {z'\y) = 5{x - y)5{x' - y'). 



(71) 



Recalling that 

{y'\z) {z'\y) 

= —F= exp iV2zj{y' -y) 



exp 



1 r 



iy-V2z'j,r + {y'-V2z'j,m72) 
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and 



G(z'; 



(Pz f d^drj 



2tt J 2tt 

we may evaluate first the integral over z/ 
dzi 



{fW)/4eii{^^n-z'^)+v{.i-z',)] ^^1^^ ^^1^,^ 



27r 



exp 



izi{V2{x-x')+T]) =5{V2{x~x')+7]) 



and the integral over z'j 



dzj 

"2^ 



exp 



iz'j{V2{y' -y)- r,)] = 5(V2 (y' - y) - v) 



Then, apart from the pre-faetor 



^e^p(-l{[y'+y''] + [x'+x''])], 



(73) 
(74) 

(75) 
(76) 



we have to evaluate the integral 

dzRdz'j^ f ^c(«'+"')/V[«(^«-^«)],5(V2(a; - x') + v)S{V2{y' - y) - v) 



2tt 

The integral over rj yields 

dr/e'''/4(5(V2 {x - x') + ij)5{V2 {y' - y) - r/) (78) 

.(ib'^'+ih'-^]') I dr^s{V2{x - x') + rj)S{V2{y' - y) - v) 



(77) 



= 6^P ( 4 y -yf + ^ W - xf j -^^(^ -x' + y' - y), 
so that the pre-faetor beeomes 

^exp(^-^{[y' + yf + [x' + xf)^^d{x-x' + y'-y), (79) 

while the gaussian integrals over z^ and z^ give 

dzi^dz^e(-2"«+^«[^(^+«')-'?])e(-24+^«[V2(x+a;')+»c]) (80) 



2 



2 exp 



\ {[y + yf + [X + xf)] ^~e/iS^x+x'-(y+y'))V-2/i 
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Collecting all the terms, the pre-factor becomes 

^5{x-y-{x'-y')) (81) 

and the last integral to evaluate is 

f ^^^i{xW-iy+y'))V2/4 ^ 2^2,5 {x + x' - {y + y')) . (82) 

J 2tt 

So, we get the expected result 

2S{x -y-{x'- y'))S {x - y + {x' - y')) = 5 {x - y) 5{x' - y'). (83) 

8.2 Checking the resolution of the unity for the photon 
number tomography 

Our aim is to check, for any allowed value of A, the matrix form of the resolution 
of unity for the photon number tomography, Eq. (j63p . which in the position 
representation reads 

^ / _f (y' \nz) {nz\x') ix\Gx (n, z) \y) = 5{x - x')6{y - y'). (84) 

With z = {v + ifi) / \2, in the position representation {y\nz) is 



dq {y^D {z) \q) {q\n) ^ j dqS {y - q - v) exp [i (pq + iJ.iy/2)] {q\n) 

= exp[i{ny- fij^/2)]{y-iy\n) , (85) 

where the n-th Hermite fmiction {q\n) is 

{q\n) - (V^2"n!)-i/2 exp(-iq2)^„(^). (ge) 

As we have shown in Ref. [12] , the matrix element of the Gram-Schmidt operator 
G\ eventually reads: 

/ ~ \ /r\ e*'^("+^/^) 

{x\G),{n,z)\y) = 4sin2 - — =^ exp [i^ (y - x)] x (87) 

\ ' ^^^ v27rismr 

exp 



\ [V ' ^^ ' ) 2sinT sinr \) 



where we have put 



^ + 1 ^r , ^ , ^ C*" + 1 . T 4 ^ . ^T 

- = e , (t A=o = Ti") ^ A = — = — icot— ; ^=4sm — . 

A- 1 'V IA u ; er _i 2 1- A^ 2 

Formulae ([M]) . ([57)1 correct the corresponding expressions given in Ref. [T^ 
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Bearing in mind the expression of {y'\nz) given by Eq.s (|85l86p and recalling 
that z = {v -\- in) /\/2, so that tPz = dpLdv/2^ we have to evaluate 

T. f^^^P My ~ ^')] kv ~ ^ \n) {n\ x' ~ v) (a:|G(^) (n, z) \y) . (88) 



Extracting the factor exp [i/i(j/ — x)\ from Lx\G\ (n, z) |j/\, integration over /x 
yields 



27r 



exp [i^i{x' - y')] exp [i^j.{y - x)] ^ 5{x' - x - {y' - y)). (89) 



It is possible to sum the series by means of the Mehler's formula, which is 
derived in the next subsection: 



Ed) 

n=0 



z\" H^{x)H„{y) 



VT 



: exp 



z^{x^ + y^) — 2zxy 



z2-l 



In our case z = e*'^, so 
z 1 



(|z| < 1, z^±l). 



1 



(90) 



-1/2 



,-ir/2 



z^ - 1 z - z 1 2 sin r ' Vl - z^ Vz"^ - z %/-2isinr 
and the Mehler's formula reads 



-ir/2 



V— 2isii 



: exp 



2 sin T 



{e- [(x' - ^)2 + iy' - ^)2] - 2(a;' - ^)(y' - i^)} 



We are left with a prefactor 

'T\ 1 e-*^/2 gJr/2 



4 sin 



• 2 /^\ -L 
m — —^ 



\/7r ^/ — 2i sin r ■\/27ri sin r ti" 



tan ■ 



(91) 
(92) 



and the integral 



di^exp — 



smr + le 



2 sin T 



{x' - ^)- + [y' - vY 



+ .(-^--)(^^--)|x 



smr 



COST 



exp 2— T 



2sinr 
The coefficient of v^ vanishes 



{x~vf + {y- vf 



.{x-v){y- v) 



smr 



(93) 



= 0. 



(94) 



Then, apart from the prefactor 



Itan^ I I COST 

-^ exp 

2 sin T 



[x" + y"]+l^[x' + y'] + 



i{x'y' - xy) 



2sinT 



smr 



, (95) 
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we have 



dz^exp ( V 
dvexp 



icosT, , ,, icosr, , i(x' + y' — (x + y)) 

{x' + y') - {x + y)- ^ ^^^ ^ 

sm r sm r sm r 



^{x +y — (a; + y))tan■ 



27^ 



tan 



:^S{x' + y - (x + y)). 



Hence we get the expected product of deha functions 

26{x' -x-{y'- y))S{x' + y' - {x + y)) - d{x' - x)5{y' - y). 
In view of this result, we may assume x' — x,y' — y in the prefactor 



exp 



[x'2 + y'^]+^[x^ + y^] + 



i{x'y' - xy) 



1, 



(96) 



(97) 



2 sinr 2 sinr 

so that we obtain just the resohition of unity of Eq. 



8.3 A deduction of the Mehler's formula in the holomor- 
phic representation 

Here we give a deduction of the Mehler's formula (1866) by using the holo- 
morphic, i.e. Bargmann-Fock, representation of the harmonic oscillator wave 
functions: 



(V'mlV'n) = 






'm\ \/n\ 
Representation of Dirac's delta function reads: 



d^i.a = 5nm, djj,^ := c "= "^ = e~ 



(98) 

,dz* dz 
2m 



f{z*) = y e^ ' /(/*)d/i., . (99) 

We use this representation of the delta function for evaluating the integral 
J{x,yX) = fe^-'-'"/^e'^''''e^y''"-'-""/^d^iJ^l,, 

(100) 



gV22;^-zV2gV2Cy^*-C'z*V2^^ 



Then we get for J(a;, y; Q the gaussian integral: 



exp 



det 



(z,z*) 



1 1 

1 e 



J 1 1 


1 


/ z 


^ 1 e 


\z* 


1 

"1^2 

exp 


i^Xv) 


J 


L 





V2{xXy) 



z 

z* 



dz*dz 



2m 



1 1 

1 C^ 



X 

Cy 



(101) 
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Eventually 



Jix,yX) 



VT^ 



: exp 



C' (x2 + y2) _ 2Cxy 



e~i 



(102) 



On the other side, remembering that the generating function of the Hermite 
polynomials in the position representation is 



V2xz-z^/2 ^ V -^" ^^^ f — \ 



E 

n=0 



\V2j ' 



(103) 



we have that Eq. (|100p also reads: 

n,m— 

Again we use the Dirac's delta and get for J(x, y; C) : 

Hn{x)Hjn{y) f „/>*xm, V^ ("' Hn (x) H^ (y) 



E 



2(n+m)/2^I^I 



z"(cz*)™d/i, = y v^ 



)/^Vr 



Ittj I 



(V'„IV'„ 



(105) 
By comparing this expression of J{x, y; C) with the previous one of Eq. p02p 
we finally obtain the Mehler's formula 



E U 



n=0 



CY Hn{x)Hn{y) 



1 



v/T^ 



: exp 



--2 t^2 



C [x' + y") - 2Cxy 
C^-1 



(106) 



The radius of convergence of the series is 1, but from the previous evaluation 
we get that the values |C| = 1, C 7^ il £^re allowed because 



det 



1 1 

1 e 



^0. 



The limits C -^ ±1 are 
1 



lim 



exp 



C' (a;2 + y2) _ 2Qxy 



e-1 



(107) 



= ^A(5 (x T y) . (108) 
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